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Global scale impacts modify the physical or thermal state of a substantial fraction of a 
target asteroid. Specific effects include accretion, family formation, reshaping, mixing and 
layering, shock and frictional heating, fragmentation, material compaction, dilatation, stripping 
of mantle and crust, and seismic degradation. Deciphering the complicated record of global 
scale impacts, in asteroids and meteorites, will lead us to understand the original planet-forming 
process and its resultant populations, and their evolution in time as collisions became faster and 
fewer. We provide a brief overview of these ideas, and an introduction to models. 


1. Introduction 

The most important parameter governing the global ex¬ 
tent of an impact is the mass ratio of the projectile to the 
target, 7 = M 2 /Mi. In the case of a cratering event this 
ratio is small, and there is a well-defined geometric locus. 
Crater scaling then becomes a powerful tool (e.g. Housen 
et al. 1983) that allows simple analytical approaches to be 
applied to determine whether an impact ‘goes global’ - for 
instance, whether the surface is shaken everywhere to the 
escape velocity, or whether the target is shattered or melted. 

At the other extreme, as M 2 —> Mi, there is no impact 
locus, so the mechanics and dynamics are complex (As¬ 
phaug 2010). Crater scaling does not apply, even though the 
impact physics is fundamentally the same. By definition, 
these similar-sized collisions (SSCs) are global events. Un¬ 
like most cratering impacts, they involve substantial down- 
range or even escaping motion of the projectile M 2 , de¬ 
pending on the impact velocity Vimp and the impact angle 

e. 

Impact velocity is the next most important parameter, be¬ 
cause of its dynamical and thermodynamical consequences. 
In order for an asteroid to be eroded in a cratering impact, 
for example, one projectile mass of material must be ejected 
to Vej > Vesc, where Vej is the ejecta velocity and Vesc is 
the escape velocity described below. The ejection velocity 
in turn scales with Vimp, so that a cratering projectile has to 
strike at a few times v^sc in a competent rocky target if it is 
to cause net escape of material, and an order of magnitude 
faster to cause net mass loss from a highly porous target 
(Housen and Holsapple 2011). This implies a very different 
impact evolution depending on porosity. 

Impact velocity also represents a specific collisional ki¬ 
netic energy Q = ^{Miv"/ + M 2 V 2 )/{Mi + M 2 ) where 


vi < V 2 are the velocities of the target and projectile in the 
center of mass frame. If the energy is sufficiently intense, 
Q > Q*g, then shattering occurs, breaking the solid bonds 
of the asteroid into pieces no larger than Mi/2. If gravita¬ 
tionally bound (ejected at < Vesc) then shattering produces 
a rubble pile as defined below; otherwise if fragments are 
escaping, the result is a collection of new asteroids. This 
is the classic example of hitting something so hard that 
you break it. Catastrophic disruption requires a collision 
of greater energy Q > Q*^, > Q*g, to break any solid bonds 
and also to overcome internal friction (see below), and to 
exceed the gravitational binding energy, thereby dispersing 
the fragments to Vej > Vesc- For massive bodies, interme¬ 
diate energy collisions Q'/) > Q > Q*g can lead to compli¬ 
cated (altered and reassembled) geologies. 

These are idealizations because impact energy is not de¬ 
posited uniformly inside a target. Much of this chapter is 
to study how this deposition occurs, and what it does. The 
impact angle 9 is especially important in this regard, es¬ 
pecially for similar colliding masses, since only a limited 
amount of angular momentum can be accreted in a colli¬ 
sion, and because objects of comparable diameter tend to 
suffer grazing collisions more often than not. And finally, 
two asteroids of masses Mi and M 2 cannot be thought of 
as colliding in isolation, even if one ignores all the other 
asteroids and planets: both bodies orbit the Sun and their 
fragments continue orbiting on intersecting orbits, so their 
interaction extends long after the original collision. 

1.1. Mass effects 

Generally speaking, impacts slower than Vesc cause ac¬ 
cretion, and impacts faster than Vesc cause erosion, with the 
specific boundary depending on the impact angle 9. This is 
effectively the case for cratering impacts and for similar- 
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sized collisions, so we consider the governing parameter 
Vimpivesc- But Vimp Can be no slower than Vesc, the im¬ 
pact velocity of two spheres falling from infinity with initial 
relative velocity Vrei = 0: 

Vesc = v/2G(Mi+M2)/(i?i+i?2) (1) 

where Ri and R 2 are the corresponding radii. As a rule of 
thumb, Vesc = Rkm, in ni/s, where Rkm is the radius of 
a spherical asteroid in km; this holds exactly true for bulk 
density p = 1.9 g/cm^. 

For a 2-body encounter in the absence of gas or other 
drag effects, conservation of energy implies that 

= \J ( 2 ) 

Concerning accretion and erosion, whether by cratering or 
SSC, it is consequently this relative velocity above Vesc that 
matters most: 

^ — Vrel!Vesc (3) 

The magnitude of (p defines the dynamical ‘kick’ that drives 

material out of the gravitational potential of the target. 

Dynamical studies of planet formation have tended to 
ignore the substantial differences between a collision (two 
objects’ radii intersecting) and an accretion. A collision at 
V = Vimp at a separation distance r < i?i -f i ?2 at contact 
angle 0 is a geophysically and astrophysically complicated 
event, so that even the slowest merger (vimp ~ Vesc, 4> ^ 0) 
is an imperfect accretion. Yet most A^-body simulations to 
date implicitly assume that whenever two planetary bodies 
touch each other, they become a single object orbiting the 
Sun. This is generally an invalid assumption (Agnor et al. 
1999). 

Accretion efficiency ^ is defined as the fraction of the 
projectile mass M 2 acquired by the target Mi. Perfect ac¬ 
cretion ^ = 1 makes a final merged body Mp = Mi -f M 2 , 
so in general 

Mi)/M2 (4) 

Accretion efficiency is always a few percent less than 1, 
even in Vrei = 0 collisions, because mergers release a sub¬ 
stantial fraction of the gravitational binding energy com¬ 
pared to two contacting spheres. Mass is further lost by 
angular momentum redistribution, flung out by the spiral¬ 
ing, merging protoplanets. Mass loss in otherwise ‘perfect 
mergers’ can lead to satellite formation (e.g. Canup 2004) 
and percentage losses of escaping remnants (Asphaug and 
Reufer2013). 

In cratering collisions, P^{(j>) is a relatively smooth func¬ 
tion, becoming negative for ~ 1, except in the case for 
highly porous bodies, where ^ may remain positive (accre¬ 
tionary) to velocities tens of Vesc in a process of compaction 
cratering (Housen and Holsapple 2011). SSCs (7 > 0.03) 
are substantially different because of geometrical effects, 
and this leads to generally low accretion efficiency. For 
SSCs accretion efficiency is a sensitive function of veloc¬ 
ity ^{(p) that starts near ^(0) < 1. At moderately off-axis 


angles (only ~ 15° — 30°) SSCs are grazing, in the sense 
that most of the colliding material does not physically in¬ 
tersect. The colliding bodies can ‘bounce’, with outcomes 
very sensitive to <p and 6. 

As reviewed by Asphaug (2010), there is an abrupt tran¬ 
sition in accretion efficiency from ^(0) ~ 1, to ^ « 0 (little 
mass change to the target, but huge transformation to the 
impactor) over the velocity range 1.2 < (p < 2 — 3. We call 
these ‘hit and run’ (see chapter by Scott et al. 2015). And 
finally, erosion and disruption (^ < 0) occur over the range 
(p > 2 — 3. The exact boundaries of these curves depend 
on the composition and degree of differentiation of the col¬ 
liding bodies, as explored in further detail by Stewart and 
Leinhardt (2012). 

One of the major unknowns of planet formation, is the 
kind of collision that dominated the accretion of mass by 
asteroid parent bodies. If it was accretion of myriad much 
smaller objects (cratering, 7 —> 0), then effective com¬ 
paction might facilitate accretion, with highly porous plan- 
etesimals acting as ‘sponges’ to sweep up much smaller ob¬ 
jects. Truly primitive bodies would contain a record of com¬ 
paction events (e.g. Belton et al. 2007). If SSCs dominated 
(7 —> 1), then small planetesimals would be unlikely to ac¬ 
crete further once random velocities were excited to faster 
than Vesc- Turbulent stirring of the nebula would lead to 
their disruption (Benz 2000), one of several ‘size barriers’ 
to primary accretion (Weidenschilling and Cuzzi 2006). 

One possible solution is to bypass these bottlenecks dur¬ 
ing the gas phase, for example by streaming instabilities 
followed by pebble accretion as described in the chapter by 
Johansson et al. (2015). Such models indicate the rapid for¬ 
mation of asteroid progenitors ~ 100 — 1000 km diameter 
in the presence of the nebula. These bodies would be well 
sorted by mass, so that their further collisional interaction 
would be dominated by SSCs. Another implication is that 
many of these bodies (those containing chondritic abun¬ 
dances of active radiogenic Al) would likely be melted by 
radionuclide decay during the timeframe of their collisional 
interaction. The surviving original bodies would include a 
preponderance of hit and run relics according to Asphaug 
and Reufer (2014), by the attrition of a population that is 
mostly accreted into growing planets. 

1.2. Thermodynamical effects 

Impact heating is also sensitively dependent on velocity. 
In addition to radiogenic heating, shock and frictional heat¬ 
ing were common throughout early solar system history, as 
recorded in meteorites, with shocks being particularly (but 
not uniquely) effective at dissipating impact kinetic energy 
Q into heat (Ahrens and O’Keefe 1994; Melosh 1989). 

Shocks are produced when the impact velocity exceeds 
the sound speed, Vimp ^ c^. Essentially, impact momen¬ 
tum is added to the medium faster than it can be transported 
away as a pressure wave, leading to a discontinuous ‘jump’ 
that causes irreversible thermodynamical effects (heating, 
compression, momentum deposition). Because impacts are 
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generally faster than Vesc, this means that planet-scale colli¬ 
sions always involve shocks (Vesc of tens of km/s), whereas 
asteroid-scale impacts must be stirred up externally to high 
random velocity to produce shocks. 

Random velocities in the Main Belt are excited by gravi¬ 
tational interactions with the terrestrial and giant planets, to 
velocities that are typically Vrei ~ 2 — 8 km/s, exceeding 
Vesc by orders of magnitude. Although we have no direct 
measurements of asteroid sound speed, an upper limit is 
probably that of the lunar mantle, Cg ~ 7 — 8 km/s. A more 
representative sound speed may be that measured by Flynn 
(2005) and others for competent fragments of chondrite me¬ 
teorites, Cg ~ 1 — 3 km/s. Thus a great deal of shock pro¬ 
cessing occurs on small bodies in the evolving and modern 
solar system. But unlike planet-scale collisions that retain 
their shock-melted materials, the shocked ejecta in asteroid- 
asteroid collisions are generally lost as interplanetary dust. 

The sound speed is much slower in rubble piles and near¬ 
surface regoliths, perhaps only of order ^ 100 m/s, so in 
this crushable zone it is possible to form melted materi¬ 
als that could be retained in place, or at least, be ejected 
at < Vesc- Housen and Holsapple (2011) find that crater 
ejection velocities in highly porous silicate material can be 
more than an order of magnitude slower than crater ejection 
velocities in competent materials. This means that in addi¬ 
tion to being more prone to shock melting and compaction 
melting, the melted material may be less prone to escaping 
on highly porous bodies. 

1.3. Velocity relationships 

Impact velocity relative to escape velocity, Uimp/fesc. 
governs accretion efficiency and its mass effects, while im¬ 
pact velocity relative to sound speed, Vimp/cg, governs the 
formation of shocks and its thermodynamical effects. Ac¬ 
cordingly there are four quadrants of global scale impacts, 
as summarized in Figure 1. 

Shock-forming impacts are to the right of the diagonal 
band, represented here by a wide possible range 0.5 km/s 
< Vimp < 8 km/s. Mass loss and catastrophic disruption 
(^ < 0) occur above the horizontal band, broadly the range 
Vimp ^ 2 to 20vesc (competent or highly porous target; 
cratering or SSC). Asteroid family formation, reviewed in 
the chapter by Michel et al. (2015), occupies the upper right 
quadrant, while parent-body accretion occupies the lower 
left. 

The transitions across quadrants (grey shaded zones) are 
just as important as the end member scenarios, perhaps 
moreso. Consider the horizontal band: most collisions will 
have Vimp equal to a few times Vesc in a self-stirred pop¬ 
ulation, because the random velocities grow to be propor¬ 
tional to the escape velocities of the largest members (e.g. 
Safronov 1972). 

As for the sound speed transition, it appears to be coin¬ 
cidence that the range of Cg for asteroid materials overlaps 
the expected collisional velocities. Random velocities are 
slower in the outer solar system, ~ 1 — 2 km/s, but so is the 


sound speed of icy bodies. So in both regions, collisions are 
expected to occur on either side of the ‘shock zone’, at least 
in the present solar system. 

During parent body accretion, most events were at or 
near this line (low velocities), and towards the left (small 
sizes). As planets grew and stirred things up, events stayed 
near the line at large sizes, but rose to high above the line 
for small sizes, grinding down that population into today’s 
generation of asteroids. 

In the first 10-100 Ma after dissipation of the dust and 
gas in the solar nebula, when collision probabilities were 
substantially higher than they are today, the average colli¬ 
sion speed in the asteroid belt grew from < 1 km/s to sev¬ 
eral km/s. Hence, impacts straddled the shock transition 
zone when impacts were most influential in defining aster¬ 
oid parent body geology. 

1.4. Overlapping timescales and linked models 

Global scale impacts involve a diverse range of overlap¬ 
ping processes and timescales, making it important to un¬ 
derstand and utilize linked models. Because asteroidal tar¬ 
gets vary from tens of meters to hundreds of kilometers di¬ 
ameter, global scale impacts can take place in less than a 
second, on the ballistic timescale 

n ~ 2r/ '^imp (5) 

where r is the radius of the projectile and Vimp is the impact 
velocity; or in seconds to minutes, on the stress propagation 
timescale 

Tg ~ 2R/c (6) 

where 2R is the target diameter and c is the sound speed (or 
shock speed); or in hours to days, on the gravity timescale 

Tg ~ yJiTlIGp (7) 

where G is the gravitational constant and p is the asteroid 
bulk density. 

In the end-member regimes (e.g. self-gravitating plan¬ 
ets without strength, or monolithic rocks without gravity) 
key variables can be eliminated and simple models or scal¬ 
ing predictions can be viable. Similarly, if one can bypass 
a complex phase or timescale of a collision using an effec¬ 
tive approximation, then a simpler analysis can be available 
- for instance, skipping the self-gravity phase by applying 
an escape-energy approximation, or replacing the hyperve¬ 
locity impact calculation with a ‘^-model’ velocity field 
derived from impact cratering experiments (c.f. Melosh 
1989). It is an open question whether these and other sim¬ 
plifications are well justified for modeling global scale as¬ 
teroid impacts, when the balance of force is so subtle. 

Another approach is to use end-to-end linked numerical 
models. For asteroids up to tens of km diameter, Tg, 
are distinct, making it is possible to consider three separable 
regimes of impact, using one calculation as input for the 
next. This of course makes the assumption that not a lot of 
physics happens in the intermediate regimes, which as we 
shall see is not always the case. 
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Fig. 1.— Global scale impacts span four quadrants in a graph of impact velocity normalized to target escape velocity, Uimp/uesc, 
versus target radius. The horizontal hand marks the approximate transition from accretionary to erosive/disruptive impacts, a few times 
Vesc- The diagonal band defines impact speeds comparable to the sound speeds for a wide range of solar-system materials. Planetesimal- 
scale accretion (hottom left) does not involve shocks, while planet-scale accretion (bottom right) does. On small bodies, mass loss and 
disruption do not require shocks. Large Vim.plv^sc. requires gravitational stirring by larger bodies, so that events to the right of the graph 
tend to be damped (no larger bodies; gravitational drag by smaller bodies) while events to the left tend to he excited, for instance the 
modem Main Belt. 
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Consider the example of forming a ~ 10 km crater on 
a ~ 20 km diameter asteroid (e.g. Stickney on the mar¬ 
tian satellite Phobos, or Shoemaker on asteroid Eros). Here 
Tb ^ 10“^ second, Ts is seconds, and Tg is hours. An im¬ 
pact hydrocode (Anderson 1987) is a versatile simulation 
tool (see chapter by Jutzi et al. 2015 and below) that can be 
run for several Tb to compute the shock pressures, energies 
and temperatures and the initial momentum distribution em¬ 
placed by the projectile. The accuracy of this step is limited 
by our understanding of asteroid target geology, which in 
even the best codes is represented simplistically (equations 
of state and crush models). For example, rate dependent 
effects are not typically included in hydrocodes. 

After the impact coupling phase is complete, shock pres¬ 
sures and stress states can be evolved using a wave prop¬ 
agation code (e.g. Richardson et al. 2004; Blitz et al. 
2009), or using the hydrocode itself within the limitations 
and caveats discussed below. Hydrocodes are not gener¬ 
ally designed for seismic wave propagation, although as 
discussed further below, the most modern codes are begin¬ 
ning to do a good job on the equation of state (EOS) as¬ 
pects of the problem, and the fracture mechanics parts of the 
problem, and also the granular physics which is intimately 
connected to the acoustics and flow of disrupted material. 
Simpler approaches are also feasible, and sometimes rec¬ 
ommended, especially those that are anchored in laboratory 
and field data. A power law of stress attenuation with dis¬ 
tance (Cooper and Sauer 1977) can be applied (e.g. Ryan 
and Melosh 1998; Asphaug 2008), as well as direct compar¬ 
ison with nuclear explosions (Asphaug and Melosh 1993). 

If Tg is orders of magnitude longer than the other 
timescales, then the earlier phases of a collision, establish¬ 
ing the pressure evolution, can be calculated with no gravity 
at all, or with gravity expressed as an overburden term and 
constant gravity. The largest scale events conclude with a 
gravitational and rotational-mechanical evolution over sev¬ 
eral Tg, as modeled by Michel et al. (2003), Durda et al. 
(2004), Leinhardt and Stewart (2012) and others; see chap¬ 
ter by Michel et al. (2015). 

While gravity can be ignored below some small size, it 
is not clear what size this is - gravity being a long-range 
force and strength, friction, and cohesion being short-range 
forces. The escape velocity of a 100 m diameter rubble pile 
is a few cm/s, so massive debris evolution on a small aster¬ 
oid can perhaps be astrophysically similar to much larger 
scale events. 

On the other hand, the central pressure inside a 100 m 
asteroid is only ~ 1 Pa (10 dyn cm“^), and while not zero, 
it is hundred times smaller than the cohesions measured for 
upper lunar regolith. That is, granular cohesion may well 
dominate over gravity, inside of asteroids up to kilometers 
diameter, as discussed further in the chapter by Scheeres et 
al. (2015). 

To make matters worse, granular materials can be cohe¬ 
sive in the absence of internal vibrations, and fluidized by 
the acoustic energy of a global scale collision, responding 
like a mobilized global landslide (e.g. Collins et al. 2004). 


Anyone who has played with sand is familiar with the solid- 
liquid behavior of granular media. Models are just barely 
catching up to our awareness of the complex physics that 
can be at play in a global scale asteroid collision. 

1.5. Small giant impacts 

For large planetary collisions a separation of timescales 
is generally not possible. In the Moon-forming giant im¬ 
pact, Tb, Ts,Tg are all ~ lO^s, so that all processes must be 
calculated simultaneously. Fortunately giant impact com¬ 
putations are greatly simplified by two facts: planets can 
be treated as idealized fluids on these timescales, and the 
computational timestep dt can be long, requiring fewer hy¬ 
drocode iterations, as described in Section 2. 

For numerical accuracy, dt must be substantially shorter 
than the wave-crossing-timescale dx/cg where dx is the 
spatial resolution and Cg is the sound speed; or in the case 
of shock-forming impacts, dx/vimp- Otherwise informa¬ 
tion travels faster than the physical propagation velocity 
in a computation. In medium resolution simulations of 
the Moon forming giant impact dt > 10s, which is large 
enough that an all-in-one 3D hydrocode approach is feasi¬ 
ble for the phases leading to the capture of protolunar ma¬ 
terials. 

Although we don’t consider giant impacts in this chap¬ 
ter, we do consider what might be called ‘small giant im¬ 
pacts’, that is, planetesimal-scale SSCs at speeds propor¬ 
tional to the bodies’ mutual escape velocity. Fike giant im¬ 
pacts, these straddle the erosion-accretion transition in Fig¬ 
ure 1. This is thought to be the regime where the parent 
bodies of asteroids were accreted (see chapter by Johans¬ 
son et al. 2015). Modeling small giant impacts is most 
challenging computationally, because of the short timestep 
combined with complex physics. Consider a 3D impact 
simulation of 10 km planetesimals colliding; dt < 0.1s 
so that it requires > lOOx more computational effort than a 
comparable Moon-forming calculation. For well-subsonic 
collisions the timestep can be sped up by artificially soft¬ 
ening the equation of state, thus reducing the sound speed 
(c.f. Asphaug et al. 2011) while not greatly changing the 
response so long as the collision remains subsonic. 

Although scale-similar to giant impacts, there is usually 
no shock when the colliding bodies are smaller than about 
1000 km diameter (Asphaug and Reufer 2014). This makes 
the physics quite different, although there is still substantial 
heating, shear alteration and melting, all ultimately com¬ 
pensating for the loss of gravitational binding energy. The 
slowest possible collisions between ~ 100 km planetesi¬ 
mals occur at ~ 100 m/s, like two high-speed trains collid¬ 
ing, so friction and compaction do a great amount of work. 
It is a cutting edge problem in computational modeling, one 
that can’t easily be remedied using linked models on sep¬ 
arate timescales, and therefore is in need of experimental 
work on representative bodies. 
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1.6. Rubble piles 

The consequences of accretionary collisions may thus 
dominate the geologic records of the largest asteroids and 
primary planetesimals. In small asteroids, on the other 
hand, this record has probably been beaten down. Be¬ 
cause their relative collisional energies are more disruptive, 
compared to their binding energies (larger Q, stirred up by 
planets), they are second generational or further, as repre¬ 
sented by asteroid family formation (chapter by Michel et 
al. 2015). 

Today most asteroid mutual collisions are erosive or dis¬ 
ruptive, Vrei ^ Vesc- Impacts produce craters, sometimes 
with hemispheric morphologies (mega-craters). If crater¬ 
ing ‘goes global’, as discussed below, the result is global 
scale resurfacing of the asteroid, by ejecta deposition and by 
shaking down of pre-existing topography. And of course, 
impacts contribute great energy of heating, melting, frag¬ 
mentation, and other forms of material transformation. 

Strength is partly a material property, and partly a prop¬ 
erty of the scale and the dynamics of the collision (Hol- 
sapple 2009). Size- and rate-dependent dynamic strength 
effects (Grady and Kipp 1980; Melosh et al. 1992; Housen 
and Holsapple 1999) make large bodies much easier to de¬ 
stroy, mechanically, than smaller ones. However, gravita¬ 
tionally large objects are harder to destroy. This means that 
it is much easier to create a rubble pile, breaking rocks to 
pieces, than to destroy it, sending half of its mass onto es¬ 
caping trajectories (reviewed in Davis et al. 2002). In the 
1990s it was demonstrated in computer simulations (e.g. 
Benz and Asphaug 1999) that asteroids larger than a few 
100m diameter are likely to be pulverized completely, sev¬ 
eral times over, before they are likely to be destroyed. But 
as noted by Davison et al. (2013), this does not mean that 
a given subcatastrophic event is inevitable - just that com¬ 
minution is much more likely than disruption, Q s ^ Qd- 

According to this analysis, asteroids larger than a few 
100 m in a collisionally evolved population become rubble 
piles. How rapidly or completely they do so, is an open 
question. It does not mean that asteroids larger than 100 m 
can be treated as liquid bodies in simulations, or that smaller 
asteroids are necessarily hard rocks (Holsapple 2009). In¬ 
deed it has become obvious in recent years that there is a 
broad transition regime, rather than a transition size, that 
spans perhaps ~ 0.1 — 100 km diameter, where mechani¬ 
cal, gravitational and impact momentum stresses need to be 
calculated simultaneously, altogether. This includes all but 
the very smallest asteroids. 

1.7. Giant craters 

Giant craters are principal probes of interior geology. 
They excavate deep inside the asteroid and produce re¬ 
verberating stresses that cause global surface modification, 
massive faulting and overturn. Asteroids with giant craters 
(most of the ~ 10 km ones evidently) must be ‘strong’ in 
some sense, in order to survive the mega-cratering events 
themselves. On Mathilde, several hemisphere-sized simple 


craters were excavated (or crushed; Housen and Holsapple 
1999) without destroying the pre-existing craters, not even 
their steep walls and rim structures (so far as we can tell at 
flyby resolution; Figure 2). For an unconsolidated material 
one might expect sustained reverberations that would trig¬ 
ger an asteroid to relax into a spheroid, but Mathilde did 
not; this large asteroid retained some of the most spectacu¬ 
lar topography in the solar system. 

To sustain this kind of dramatic topography, one might 
expect that a solid target structure is required. But this 
presents a paradox: breaking holes this large into a com¬ 
petent asteroid, will catastrophically disrupt the asteroid, 
according to strength-scaling rules and hydrocode models 
with strength (Melosh et al. 1992). So we get backed into 
a corner, that these large craters cannot have formed in the 
strength regime. This does not automatically mean they are 
gravity regime craters, as we shall see. 

Of known asteroid craters, only Rheasilivia and Ven- 
nenia, the two largest on Vesta (Schenk et al. 2012; Russel 
et al. 2015, this volume) show clear evidence for gravita¬ 
tional rebound and central peak formation. These two giant 
basins overlap in Vesta’s southern hemisphere, the ~ 500 
km diameter Rheasilvia basin overprinting the earlier ~ 400 
km Veneneia, offset by some 40°, on an asteroid only some¬ 
what larger (~ 530 km average diameter). Like Mathilde, 
Vesta provides an example of the interaction of younger and 
older giant craters overlapping, but in this case their inter¬ 
action is notable. Unlike Mathilde it is easy to recognize 
the sequence of events, in part because of the much more 
detailed global imaging but also due to the fundamentally 
different cratering mechanics on a body 10 times larger. 

A number of past (Asphaug 1997) and recent (Jutzi and 
Asphaug 2011; Ivanov and Melosh 2013; Bowling et al. 
2013; Stickle et al., 2015) numerical impact modeling stud¬ 
ies of Vesta have investigated the formation of its mega¬ 
craters and their resulting global effects, and the potential 
for producing ejecta (the V-type asteroids; Binzel and Xu 
1993). Modeling two overlapping mega-craters is challeng¬ 
ing, and scientifically interesting, because of Vesta’s non- 
spherical shape and gravity potential, rapid rotation, and 
differentiated structure. In addition, a 3D approach is re¬ 
quired, since the problem is strongly non-axisymmetric. 

Jutzi et al. (2013) studied the sequential formation of 
Veneneia followed by Rheasilvia using a 3D Smooth Par¬ 
ticle Hydrodynamics (SPH) impact code described below. 
The goal of these simulations was to start with a non¬ 
rotating differentiated d = 550 km sphere, impact it once 
(to form Veneneia), then place the resulting mega-crater on 
axis, spin it to period Pvesta ~ 5.3 h, and finally impact 
the rotating cratered body 40° off-axis to end up with the 
observed topography of Vesta. Thus, Vesta’s topography 
is created by a complex process of ejection and deposition, 
requiring a realistic treatment of the ejecta gravitational dy¬ 
namics and the post-impact rheological response. A pres¬ 
sure dependent strength model, a tensile fracture model, a 
friction model, and self-gravity were included. Finally, to 
‘freeze in’ the central peak structures, a block-model ap- 
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Fig. 2.— NASA’s NEAR Shoemaker mission obtained flyby images of the 50 km diameter asteroid 253 Mathilde (Veverka et al. 
2000). The image at left was obtained as the spacecraft approached in the direction of the Sun; the image on the right is from greater 
distance at much lower phase angle (43°). Mathilde is a great puzzle for the study of global scale impacts: each giant crater should have 
reset the asteroid’s global shape, according to simulations. Instead each crater seems to have formed in isolation, without degrading the 
pre-existing craters. If coming in at random angles, these successive global scale impacts would have spun Mathilde into fast rotation, if 
even a fraction of the incoming momentum was accreted (Asphaug et al. 2002). But the asteroid is one of the slowest rotating bodies in 
the solar system, P^ot =418 hrs. 


proximation of acoustic fluidization was used. The result 
(Plate ??) is in good agreement with the shape and topogra¬ 
phy of Vesta (Russell et al. 2015, this volume). 

Understanding the formation of Vesta’s giant craters and 
the provenance and specific distribution of ejecta is key to 
understanding the observed properties of Vesta, such as the 
topography and the surface mineralogy. The proper in¬ 
terpretation of surface materials (and exhumed V-type as¬ 
teroids and HED-type meteorites) lead to important con¬ 
straints on models of internal structure (McSween et al. 
2013; Clenet et al. 2014). Jutzi et al. (2013) found that 
a significant fraction of the rocks exposed in the South Pole 
region should come from >50 km deep, and up to ~60- 
100 km deep in the central mound of Rheasilvia and in 
the region where the two basins overlap. Also, according 
to the model, a significant fraction of the material that es¬ 
capes Vesta during the more recent impact, which would be 
the dominant source for HEDs, comes from greater than 
40 km deep (Clenet et al. 2014). This is puzzling be¬ 
cause olivine-rich mantle rock has not been detected in the 
Veneneia/Rheasilvia region on Vesta, and there is also a 
lack of mantle samples among the HED meteorites. 

The rapid rotation of Vesta and other asteroids can lead 
to curious mega-cratering results. In the course of reproduc¬ 
ing Rheasilvia on an already-spinning Vesta, Jutzi and As¬ 
phaug (2011) found that ejecta from mega-craters would be 
draped back on the asteroid in complex overlapping lobes, 
falling back in the rotating frame. Since for Vesta the crater 
rebound timescale is comparable to the rotation timescale. 


the Rossby number of the flow i?o ~ 1, and the result is ac¬ 
tually ‘Coriolis topography’ as shown in Eigure 3 - that is, 
crater rebound morphologies that form a spiral pattern rel¬ 
ative to the impact locus. This provides a potential window 
to further explore mega-cratering mechanics, especially in 
materials not covered beneath impact melt (Schenk et al. 
2012). 

1.8. Seismic erasure 

There is a strong correlation between the normalized di¬ 
ameter of the largest undegraded crater D^ax seen on an 
asteroid, and the diameter D of the asteroid. Dmax is well 
defined for most imaged asteroids, and for a number of the 
best resolved radar-imaged asteroids (see chapter by Ben¬ 
ner et al. 2015). Define x = Dmax/D. The smallest as¬ 
teroid seen up close in high detail, Itokawa (mean diameter 
D = 330m) has a largest crater Dmax ~ 30 m, that is, 
X ~ 0.1. The crater is not obvious to see. Eor Mathilde 
(D ^ 53km), X ^ 0.6 - and is reasonably well determined 
even from medium-resolution images. Eor Eros (D ^ 20 
km). X ~ 0-4- The trend of x varying inversely with D 
appears to hold generally (Asphaug 2008). Either global 
spanning craters do not form on small asteroids, or they are 
unable to retain them. 

Shoemaker Crater on Eros shows a different kind of ev¬ 
idence (and so perhaps would Mathilde if imaged at good 
resolution). The formation of this 7 km recent basin appears 
to have shaken down regional pre-existing ~ 100 m scale 
topography, both next to and opposite the crater (Thomas 
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Fig. 4.— Simulating the opening and collapse of the Rheasilvia complex crater on asteroid Vesta; from Jutzi et al. (2013). The rapid 
rotation of the asteroid leads to strong Coriolis forces in the rebounding ejecta, as indicated by the vectors and in agreement with the 
patterns on the crater floor. Outward vectors plot the distant and escaping ejecta. 


and Robinson 2005). The seismic erasure may have been 
impulsive (as seems to be the evidence), taking seconds to 
minutes as the initial pulse radiated from the impact; or else 
diffusive as acoustic energy remained trapped inside the as¬ 
teroid for many Ts . 

If seismic erasure is efficient, then why has not whole¬ 
sale seismic degradation shaken down the early giant craters 
of Mathilde? One approach (Asphaug 2008) is to define a 
self-consistent critical crater diameter Dcrit for which the 
above is true, and to show that Dcrit increases with asteroid 
diameter D. Craters larger than Dcrit shake down all pre¬ 
existing topography smaller than Dcrit- The giant crater 
that is forming, continues its excavation long after the seis¬ 
mic impulse attenuates, at least in this scenario, because 
the crater formation timescale in the gravity regime is many 
times Ts- A crater forms as a solitary feature if it has shaken 
down pre-existing craters of its size. Craters smaller than 
Dcrit do not cause global resurfacing at their own scale. 
Craters larger than Dcrit stand out as solitary undegraded 
craters. In this framework one can interpret that Dcrit has 
never been achieved on Mathilde (there is no one solitary 
crater), and that Dcrit would require something larger than 
Shoemaker Crater to globally erase the craters of that scale 
on Eros. 

A normalized critical crater diameter Xcrit = Dcrit/D 
can be solved as a function of asteroid diameter and tar¬ 
get properties (gravity, density, and scaling constants; c.f. 
Housen et al. 1983). Assume that the seismic impulse cre¬ 
ated by an impacting projectile decays in amplitude with 


distance r from the impact. This impulse is required to 
shake down topography to a scale equal to the diameter of 
the crater that is forming, Dcrit, out to a maximum distance 
r — D. If seismic decay with distance takes the form of a 
simple power law (e.g. Cooper and Sauer 1977), then par¬ 
ticle velocity Vp oc r~^. Crater diameter Dcrit is then cal¬ 
culated (Asphaug 2008) according to gravity regime crater 
scaling, and is substituted into the above, thus solving for 

Xcrit - 

For physically reasonable target parameters, Xcrit is pro¬ 
portional to D, accounting for the trend in crater sizes. 
On a case by case basis, matching the largest undegraded 
craters observed on spacecraft or radar-imaged asteroids 
to their theoretical maximum, does a good job at fitting 
the data if one assumes a seismic attenuation with distance 
P ~ 1.2 — 1.3 for most asteroids. While a power-law atten¬ 
uation is only an approximation to seismic dissipation, this 
value is comparable to seismic attenuations of ~0.1-10m/s 
particle velocities in unconsolidated terrestrial materials. 

For Mathilde-sized asteroids, the predicted Dcrit is 
larger than the asteroid, consistent with the handful of 
non-interacting craters greater than half the asteroid diam¬ 
eter. For Itokawa, assuming the same seismic attenuation 
as Mathilde, it turns out that Xcrit = 0.1 - the small crater 
seen on Itokawa would have resurfaced the asteroid down to 
that scale. It is an inexpensive if approximate kind of seis¬ 
mology: an Itokawa-sized asteroid exhibing larger craters, 
would have substantially greater seismic attenuation. The 
different behavior is in this case not because of material 




Fig. 3.— a, Asteroid Vesta as seen by Dawn; b, final result of the 
simulation of two large impacts in the southern hemisphere; and c, 
Lambert azimuthal projection (equal area) of the southern hemi¬ 
sphere in the model. Shading indicates the elevation (in km) with 
respect to a reference ellipsoid; for color see Jutzi et al. (2013). 
Photo credit for a, NASA/JPL/UCLA/MPS/DLR/IDA. 

differences, but because of the differences in scale: Itokawa 
next to Mathilde, is like a refrigerator next to the Empire 
State Building. 


2. Numerical Modeling of Global Impacts 

When extrapolations lead to curious results, one turns to 
direct experimentation; when that is unavailable one resorts 
to the brute force of computer simulations. These tech¬ 
niques are featured in the chapter by Jutzi et al. (2015), 
and their application to global scale impacts is summarized 
here. These techniques include hydrocodes (reviewed in 
Anderson et al. 1987), particle codes (e.g. Richardson et 
al. 2002), and linked models. These codes rely on models, 
implemented as part of a time-stepped evolution, describ¬ 
ing the material behavior - such as friction, compaction, 
viscoelastic properties, and EOS. 

Hydrocodes integrate the hydrodynamic equations de¬ 
scribing the conservation of mass, momentum, and energy. 
They relate intrinsic scalar variables - the density p, internal 
energy U and pressure P, at position x,y,z - to the veloc¬ 
ity field Vx,Vy, Vz as determined by accelerations (forces). 
They integrate a time-evolving solution to a set of partial 
differential equations. If deviatoric stresses are considered, 
then body forces are evolved according to the stress ten¬ 
sor, introducing six additional variables in 3D; in practice a 
simplified approach is often made. 

The resulting partial differential equations are integrated 
forward in time (Benz 1990, Collins et al. 2013) according 
to accuracy requirements on the spatial resolution and the 
timestep dt. Time-evolution is cast in either a Lagrangian 
reference frame, that follows the material, or in an Eule- 
rian frame that is fixed in space, advecting material prop¬ 
erties (density, momentum, energy, composition) through a 
grid. These have advantages and disadvantages; for more 
detailed background and benchmark comparison see Pier- 
azzo et al. (2008). 

2.1. Smooth particle hydrodynamics 

In smooth particle hydrodynamics (SPH) the continuous 
fluid is represented as a Lagrangian set of particles (Benz 
1990; Monaghan 2012) that move with the flow. The prop¬ 
erties at particle i are averaged over its ~ 50 nearest neigh¬ 
bors, weighted according to their distance from particle i 
relative to the smoothing length h. So h is the grid resolu¬ 
tion. Particle density, for example, is the weighted contri¬ 
bution of mass over its neighbors (hence ‘smoothed’). The 
mass and momentum conservation equations are approxi¬ 
mated to 2nd order as kernel sums. The most commonly 
used SPH code in asteroid research is by Benz and As- 
phaug (1994, 1995); this includes elastic strength and frac¬ 
ture damage in addition to the hydro-equations. This code 
has been modified by Jutzi et al. (2008, 2013) to include 
friction and compaction. 

The SPH method is mostly applied in 3D, requiring 
million-particle calculations to resolve colliding bodies ~ 
100 particles across. A fast parallel cluster is needed for 
production investigations. Nonetheless SPH is popular be¬ 
cause it is relatively simple and numerically robust. It is 
straightforward to add new physics, and the purely La¬ 
grangian description (all variables co-located at particle 
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centers, co-moving with the fluid) makes it easy to hand- 
off to other codes, visualizations, and analytical reductions. 
Unlike the grid-based codes described below, the problem 
domain is completely described by the particles themselves, 
so there is no need to grid up empty space into which the 
problem can evolve. 

2.2. Grid-based codes 

For computation of shocks, and for resolving bound¬ 
ary conditions and abrupt material contrasts, a grid-based 
calculation generally has greater fidelity than a particle 
mesh. However, in grid-based Lagrangian hydrocodes 
mesh-entanglement can result. Eulerian hydrocodes avoid 
mesh entanglement by maintaining a fixed grid, optimized 
for the problem. Each time step advects density, internal 
energy, momentum, and material type from cell to cell, ac¬ 
cording to the velocity and pressure differential. Material 
histories and interface positions are recorded using ‘tracer 
particles’ whose position moves with the mass flux from 
cell to cell. 

Advection of material type is problematic. Imagine iron 
flowing into a pure olivine cell, whose composition has just 
become 1% iron. What EOS do we use in the next timestep? 
It is not trivial: for instance what if one material is vapor at 
the computed pressure and temperature, and the other is liq¬ 
uid? There are many approaches to handling mixed materi¬ 
als in a cell. As for the mass flux from a mixed cell, do we 
assume 1% iron applies instantaneously everywhere in the 
cell (faster than any physical velocity), or do we concentrate 
iron near the donor cell, holding it ‘upstream’ according to 
some criterion? Eulerian codes usually require techniques 
for limiting this artificial diffusion. 

Grid-based codes have one overwhelming advantage 
compared to SPH: they allow accurate calculations in 2D 
(axisymmetric or planar) or ID (planar or spherical) sym¬ 
metrical geometries. High definition 2D computations, for 
instance where a spherical asteroid is impacted on axis, can 
be done on a modest workstation, enabling careful research 
exploring many kinds of hypotheses. Parameter exploration 
in 3D is much more limited. 

One popular grid-based code for asteroid collisional 
modeling is iSALE (Collins et ak, 2004; Wiinnemann et 
al., 2006), successor to SALE (Amsden et al. 1980, Melosh 
et al. 1992, Ivanov et al. 1997). This code is maintained 
in a community forum, and includes both 2D and 3D ver¬ 
sions and associated analytical and visualization tools, and 
several of the most relevant EOS models. Another popular 
code is CTH (McGlaun et al. 1990, Crawford et al. 2006), 
although like other sophisticated defense-related packages 
its foreign access restriction can limit its utility in academic 
collaborations. 

Which hydrocode to use depends on the problem. An 
intrinsically axisymmetric problem, such as impact crater¬ 
ing within about 45° of vertical, can be prototyped in a 2D 
code before serious computer effort is expended looking at 
the 3D details. If projectile-target contact and compression 


are of greatest interest, then a 2D or even ID planar calcula¬ 
tion will allow for the highest resolution. If the downrange 
momentum of massive ejecta is of key interest, then an ax¬ 
isymmetric model at any resolution is almost useless. 

To study the long term self-gravitational interaction fol¬ 
lowing mega-cratering or disruption, the hydrocode simula¬ 
tion can be handed off to a particle method such as we now 
describe. 

2.3. Discrete element models 

At low velocity, impacts can be modeled approximately 
according to Newton’s laws using discrete-element methods 
(DEMs). These do not model shocks or stress waves, but 
represent collisions by piles of immutable, self-gravitating, 
idealized solids with contact and restitutive forces (e.g. 
PKDGRAV; Stadel 2001; Richardson et al. 2002). Soft- 
spheres codes (Sanchez and Sheeres 2011; see chapter by 
Murdoch et al. 2015) further represent contact forces ac¬ 
cording to a viscoelastic ’dashpot’ model (not to be con¬ 
fused with SPH). 

Resolution in a DEM does not mean the same thing as 
it does in a hydrocode. In a hydrocode, the finer the res¬ 
olution, the more accurate the solution. In a DEM there 
is a correct resolution for the problem, because every ele¬ 
ment represents a discrete body. Infinitesimal particle size 
in a DEM results in infinite strength, and for physical rea¬ 
sons: fine powders are strongly cohesive due to their high 
specific contact area. In practice, however, DEMs struggle 
to resolve the physical scale; a 10® element simulation of 
Itokawa would have a finest block size of several meters. 
While sufficient to mimic the granular nature of a rubble 
pile, the contact area per unit volume is a billion times less 
than if we had a finest block size of centimeters, meaning 
that we must represent the surface physics parametrically. 

Polyhedral rubble piles codes (Korycansky and Asphaug 
2006; Movshovitz et al. 2012) are more challenging to im¬ 
plement than spheres-based DEMs due to the complexity 
of contact detection and the variety of contacts (face, ver¬ 
tex, edge), and have even worse resolution. But the phys¬ 
ical insights are important. Polyhedral DEMs exhibit sub¬ 
stantially different behavior than sphere-based DEMs be¬ 
cause of the dilatancy and grain locking effects that are only 
partially captured by spheres-based approaches. Eor exam¬ 
ple, Korycansky and Asphaug (2009) find a systematically 
higher value of Qij in planetesimal collisions at tens of m/s, 
compared with results derived using spheres-based DEMs 
(Leinhardt and Stewart 2009). 

Eor collisions faster than ~10 m/s, DEMs are not ap¬ 
propriate without special handling, because they do not ac¬ 
count for realistic mechanical behavior and thermodynam¬ 
ics effects. Comparable to a train wreck, or the collapse 
of a skyscraper, large planetesimal collisions are transfor¬ 
mative events involving fragmentation, crushing, shear fric¬ 
tion, and shocks - they are not piles of bouncing marbles. 

Although the science is still in its early stages, a linked 
model (e.g. Michel et al. 2003; Leinhardt and Stew- 
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art 2012) has proven versatile, where SPH (or another hy¬ 
drocode) evolves the high velocity collision for a few Tg, 
which is handed off as particle positions and velocities to 
a DEM. But even if the hand-off is perfectly achieved, we 
caution, as shown below, that the result may depend on un¬ 
known aspects of asteroid dynamical rheology. 

2.4. Strength and damage 

Strength describes a material’s response to volume- 
conserving deformation and to tensile pressure. Volumetric 
compaction is included in the EOS and the porosity model, 
described below. The strength model can be simple, requir¬ 
ing no material history, or can be a complicated evolution 
of strain, strain rate, temperature and pressure. 

Asteroid strength modeling used to be based around con¬ 
cepts of brittle fragmentation and fracture damage. While 
these approaches are still valuable, asteroid strength mod¬ 
eling now revolves more around concepts of cohesion, fric¬ 
tion and compaction - the subtle strengths of granular rub¬ 
ble piles and regoliths, in microgravity. 

A small strength or cohesion can exert a major influence 
on the asteroid impact process, especially during the final 
stages that govern how the feature will appear geologically, 
in spacecraft images. Cohesion as small as that of fresh 
dry snow can dominate global scale events on ~ 100 m 
asteroids, making the difference between a cratering event 
‘going global’ or remaining local. 

A shear strength model dependent on confining pressure 
P is employed in iSALE (Collins et al. 2004) and in SPH 
(Jutzi et al. 2013, Jutzi, 2014). In this simple approach the 
shear stress is limited by a yield stress, which increases with 
P due to grain interlocking and friction; 

= ^0 + 1 ^ 

1 + HiP/iYM — Yq) 

where Yq is the shear strength at P = 0 (measured in the 
lab), Ym is the asymptotic shear strength at infinite pres¬ 
sure (the Hugoniot elastic limit), and is the coefficient of 
internal friction. 

In elastic solid materials, a damage D can be accumu¬ 
lated either as a tensor Dij or a scalar (isotropic damage). 
Damage reduces the deviatoric stress, leaving behind (when 
D = 1) only the granular cohesive properties. To account 
for thermal softening, Y may be further reduced by a factor 
1 — U/Umeit where Umeit the specific melting energy, so 
that a molten material has zero strength. 

An approach used in iSALE is to accumulate shear dam¬ 
age as a non-linear function of the equivalent plastic strain 
(Collins et al., 2004). There is no rate dependence to this 
function; rate-dependent fragmentation requires a dynamic 
fracture model (Melosh et al. 1992; Benz and Asphaug 
1994, 1995). 

Completely damaged material has some cohesion; lunar 
regolith has Yb ~ 10^ P^- Cohesion of only tens of Pa 
can dominate the structural geology of small asteroids (As¬ 
phaug 2009; Holsapple 2009; Sanchez and Scheeres 2011). 


This cohesion can be further weakened by granular dynam¬ 
ical processes as discussed in the chapter by Jutzi et al. 
(2015) and by Kenkmann et al. (2013). 

Acoustic fluidization has been simulated using a block 
model (Collins 2004) to freeze in the central peaks of the 
largest craters on Vesta (Jutzi et al. 2013). Although it has 
not yet been implemented so far for small asteroid colli¬ 
sions, the attenuation of random motions responsible for 
fluidization might lock in global-scale shapes and structures 
of rubble pile asteroids. 

2.5. Size- and rate-dependence 

The number of flaws n per unit volume V in a solid is 
often described as a power law, 

n(e) = fee”" (9) 

In this equation a few flaws n will become active in V for 
a given volumetric strain e. There will be one single weak¬ 
est flaw, emin = so that an asteroid’s strength de¬ 
creases with size, as the power (Earinella et al. 1982) 

for m ~ 6 (a typical value). This is the static threshold; but 
dynamic strength increases with strain rate in a closely re¬ 
lated manner. The overall effect is a negative slope when 
plotting Q*g versus R, where R is asteroid radius - larger 
asteroids are weaker. 

On the other hand, larger asteroids are harder to disrupt, 
because the gravitational binding energy per unit volume 
(equivalently, the internal hydrostatic pressure ^ Gp^R^) 
is proportional to the escape velocity squared. Setting grav¬ 
itational binding energy equal to the rock strength Y gives 
a transition diameter above which it is more likely for an 
asteroid to be beaten into rubble many times over, than to 
be dispersed. 

One of the remarkable discoveries of the previous 
decade, was that the strength-gravity transition would oc¬ 
cur around 100-300 m diameter asteroids. This result was 
predicted in simple analyses and is confirmed in advanced 
numerical models, but the implications are still the sub¬ 
ject of considerable debate, as reviewed in the chapter by 
Michel et al. (2015). 

2.6. Equations of state 

The hydro-equations are closed using an equation of 
state (EOS) relating energy, density and pressure. Eor im¬ 
pact modeling, the EOS defines the material response to 
changes in volume, often in collaboration with a porosity 
model. The EOS is called by the hydrodynamics integrator 
(sometimes iteratively) as P{p, U) to obtain the ‘right hand 
side’ terms that drive evolution of the momentum equation. 

One popular model is ANEOS (Thompson and Lauson 
1972; Melosh 2007), whose thermodynamics variables are 
derived from the Helmholtz free energy. As a result, all 
data in an ANEOS table are thermodynamically consistent, 
something that is especially important in the detailed mod¬ 
eling of shocks. ANEOS returns temperature, pressure, and 
state of the material, mapping polymorphic and liquid/solid 
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phase transitions, specific heats, and vapor/liquid fraction. 
But thermodynamical consistency is not the same as accu¬ 
racy, and Kraus et al. (2012) find that the shock vaporiza¬ 
tion energy of silicate rocks may be a factor of 2-3 lower 
than reported by ANEOS, of fundamental importance to the 
study of high velocity events and giant impacts. 

Another popular EOS model is that of Tillotson (1962), 
effectively a stitched-together polytrope defining a cold 
condensed state and a hot expanded state. With only five 
parameters Tillotson is easily modified to new materials; 
one avoids the tabular precomputations and interpolations 
and possible version-inconsistencies that can plague calcu¬ 
lations. Tillotson is poorly suited to problems where vapor- 
melt interactions are expected to dominate, such as the post¬ 
impact protolunar ‘flying magma ocean’ (Stevenson 1989). 
Even so, Canup (2004) finds that Tillotson gives compara¬ 
ble results to more sophisticated EOS models for the initial 
characteristics of the protolunar disk. 

There is a practical advantage to Tillotson: it is simply 
defined, easily computed, and easily understood in a calcu¬ 
lation. Because it derives from the particle-velocity shock- 
velocity relation (Melosh 1989) it is highly suitable to pro¬ 
ducing the post-shock particle velocities, and that is what 
is most important in most asteroid collisions. Melt is not 
produced in Tillotson (there are no ‘phases’), but melt and 
vapor produced in asteroid collisions is usually ejected to 
escaping velocity anyway, as discussed further below. 

2.7. Porosity 

The realization that asteroids and planetesimals contain 
a large volume of pore space (e.g. Housen and Holsapple 
1999; Consolmagno and Britt 2002) has driven a recent ad¬ 
vance in modeling dissipative properties of porous media 
(e.g., Wiinnemann et al., 2006; Jutzi et al. 2008; Collins 
et al., 2011). To the extent that resolution allows, macro¬ 
porosity can be modeled explicitly using solid materials de¬ 
scribed above, as cavernous voids (Asphaug et al. 1998; 
Jutzi et al. 2008). But the idea of a loosely-packed dust rich 
structure for asteroids (Asphaug 2009) requires porosity to 
be modeled at the sub-resolution, in addition to the com¬ 
pressibility of a fully-compacted material that is defined by 
the EOS. 

Porosity is associated with a crushing strength. One ap¬ 
proach (the P-alpha model; Herrmann, 1969; Carroll and 
Holt, 1972) computes a distention a = psjp, where p is the 
bulk density of the porous material and ps is the density of 
the non-porous matrix corresponding to the EOS. The pres¬ 
sure P and internal energy U in the porous material, at the 
grid resolution, is what is needed by the hydrocode, but the 
EOS only knows how to compute these for the non-porous 
material, Pg and internal energy Ug. 

The P-alpha model uses the approximation 

P = -Pg(pg, Ug) = -Pgiap, U). (10) 

a a 

In this manner the EOS model for non-porous rock is con¬ 
verted into an EOS model for porous rock, whose pressure 


gradients at the grid scale then drive material accelerations 
and shocks. The distention varies with pressure according 
to a ‘crush curve’ derived from uniaxial laboratory experi¬ 
ments, as described in the chapter by Jutzi et al. (2015). 

2.8. Self-gravity 

Eor modeling global scale impacts into small asteroids, 
self-gravity can be neglected during the impact phase and 
applied later, since Tg ^ Tg. Otherwise, self-gravity can 
be approximated using a constant radial gravity field inside 
a spherically symmetric asteroid. The modern hydrocodes 
described above can also calculate gravity explicitly (e.g. 
Jutzi et al. 2013; Davison et al., 2013), although this adds 
considerably to the computational expense. 

Pre-compression must be included along with self¬ 
gravity, otherwise the first phenomenon is the asteroid 
shrinking and overshooting its central pressure. Eor small 
asteroids an analytical pre-compression is acceptable (As¬ 
phaug and Melosh 1993). Eor larger asteroids or more 
sensitive studies, gravitational stresses must be iteratively 
solved prior to the impact simulation. 

Eollowing a global scale collision, gravitational evolu¬ 
tion must be computed over the orbital timescale of days. 
This is many thousands of hydro-timescales at the resolu¬ 
tion of an asteroid, and thus prohibitive. This evolution 
is best computed after pressures and velocities have equi¬ 
librated, by mapping the hydrocode results into N-body 
codes or self-gravitating OEMs as described above. 

3. Thermal effects 

Eollowing its accretion in the first few Ma, a typical 
~200 km diameter meteorite parent body likely experi¬ 
enced several hundred impacts by asteroids larger than 300 
m during the first 100 Ma of solar system history. This is 
about half the impacts expected during its lifetime (Davison 
et al., 2013). Over the same period, impact speeds increased 
from an average < 1 km/s to as fast as ~10 km/s (O’Brien 
and Greenberg 2005). 

According to Davison et al. (2013), post-accretion col¬ 
lisions were only capable of catastrophically destroying at 
most a few percent of the finished asteroids. Of the remain¬ 
der, most were subjected to > 1 collision with an object 
1/20 its diameter, and a few percent experienced a collision 
with an object 1/5 its diameter. 

It is therefore no surprise that evidence for impact pro¬ 
cessing is common in meteorites (e.g. Sharp and DeCarli, 
2006), including shock effects, compaction, and frictional 
heating. However, as we now review, impact heating ap¬ 
pears to be inefficient at the global parent body scale, be¬ 
cause the same violent impact that heats and melts the ma¬ 
terial, ultimately causes the asteroid to disperse (Q > Q*jg). 
The fate of dispersed shocked material is another question 
- whether it finds its home back to the target asteroid, or 
another asteroid, or becomes ground into fine dust and dis¬ 
appears, or is sampled as meteorites. 
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3.1. Global impact heating of asteroids 

Widespread heating and melting of asteroids is evident in 
meteorites (e.g. Bizzarro et al. 2005) and one of the primary 
debates has been whether the source of this heating is by 
early accretion of radiogenic Al, or by impact processes. 

Keil et al. (1997) combined simple theoretical arguments 
with observational evidence and numerical impact simula¬ 
tions to conclude that the globally averaged effect of im¬ 
pact heating on a non-disrupted asteroid was equivalent to 
a global temperature increase of less than a few degrees. 
Others (e.g., Rubin 1995) maintained that shocks from high 
velocity collisions, especially into porous primitive targets, 
could be the source of melting, and that the meteorite col¬ 
lection itself would be biased to include the more highly 
shocked and melted surface materials. 

Davison et al. (2013) revisited the question of impact 
heating, using a Monte-Carlo simulation that accounted for 
the effect of asteroid porosity and impact velocity on the 
disruption threshold (Jutzi et al., 2010; Housen and Holsap- 
ple, 1990) and for the evolving impactor size- and velocity- 
frequency distribution (using a simulation by O’Brien and 
Greenberg 2005). For a 200-km diameter, 20% porous as¬ 
teroid, Davison et al. (2013) found that the maximum spe¬ 
cific impact energy that could be imparted, without disrupt¬ 
ing the target, is ~0.1-0.7 MJ/kg. This corresponds to a 
globally-averaged temperature increase of 100-700 K, de¬ 
pending on the impact velocity, substantially greater than 
the Keil et al. (1997) estimate. 

However, an impact just below the catastrophic thresh¬ 
old is not representative of the largest collision that is ex¬ 
pected. This just gives the maximum heating that is theo¬ 
retically possible, without also dispersing the heated mass. 
Davison et al. (2013) followed the collisional evolution dur¬ 
ing the first 100 Ma of solar system formation, assuming 
an approximately inverse square size-frequency distribution 
of impactors. They found that in most cases, the most en¬ 
ergetic impact experienced by a 200-km diameter asteroid 
that was not disrupted, was only a few percent of Q*jy. They 
concluded that significant impact heating on a ~200-km di¬ 
ameter parent body is unlikely, and in the end obtained a 
result similar to Keil et al. (1997), that in the most typi¬ 
cal scenario the maximum specific impact energy delivered 
to a non-disrupted asteroid amounts to a globally-averaged 
temperature increase of only a few degrees. 

Larger asteroids have much greater specific binding en¬ 
ergy and can therefore withstand collisions that deposit far 
more heat. An analogous situation is true of planetary cra¬ 
tering, where the largest craters are the most melt-rich. But 
a maximal event is less likely to occur on a larger target than 
on a smaller target, for an approximately inverse square size 
distribution. On this basis Davison et al. (2013) found that 
both the maximum specific impact energy and the total im¬ 
pact energy were rather constant for asteroids 100-500 km 
in diameter, a few tens of degrees of heating at most. These 
estimates assume that all of the heated materials of impact 
are retained in the non disrupted asteroid. A fraction of the 


hottest material actually escapes the asteroid, so that these 
global heating estimates are upper limits. 

3.2. Localized and buried heat 

These globally-averaged estimates suggest that short¬ 
lived radionuclide decay (^®A1) was a far more effective 
global heat source than impact. However, impact heating 
is transient, local, and non-uniform. In events where the 
projectile is much smaller than the target, impact energy is 
converted into heat within the proximal, near-surface vol¬ 
ume that is highly shocked, a few projectile radii, and zero 
appreciable heating occurs far from the impact. For small 
asteroids with low escape velocity, much of this melted ma¬ 
terial escapes, while on large asteroids it can be retained, 
requiring a different sort of calculation. 

Temperature increases attenuate with shock pressure, 
rapidly with distance, and so do ejection velocities. Only 
in rare, disruptive collisions between objects of similar 
mass can impact heating be considered to be nearly uni¬ 
form (Davison et al., 2010). The accretionary collisions 
described above (‘small giant impacts’) are global in con¬ 
sequence, but too slow to result in shock heating. Frictional 
and compaction heating during early accretion (the first ~ 3 
Ma) may be of global consequence, but are thought to be 
insignificant in magnitude compared to the prevalence of 
26 Al. 

Using the same dynamical and population framework 
described above, Davison et al. (2013) computed the vol¬ 
umetric significance of impact heating on 200-km diame¬ 
ter meteorite parent bodies during the first 100 Ma of Solar 
System history (Figure 4). Less than 1% by mass of a typ¬ 
ical undisrupted parent body experiences heating of more 
than a few hundred K, volumetrically equivalent to a sur¬ 
face layer < 300-m thick. However, of the few percent of 
asteroids that are disrupted, most experience a near-global 
temperature increase of a few hundred K, and one fifth are 
expected to undergo (near) complete melting. 

The main influence of asteroid porosity on impact heat¬ 
ing is to modify the distribution of heat. Shock heating is 
more extreme, but much more localized, in porous body 
collisions compared to impact on a nonporous body. More 
shock energy is absorbed close to the impact site on a 
porous asteroid by permanent crushing of pore space; the 
excess heating from pore collapse leads to higher proximal 
temperatures (Asphaug et al. 2002; Davison et al., 2010). 

The localized heating by sub-catastrophic collisions, and 
the distribution of heated materials beneath and around the 
crater, imply that a single impact can produce materials with 
a remarkable range in peak temperature and cooling rates. 
According to Davison et al. (2012), impacts on porous as¬ 
teroids can produce large amounts of melting at speeds as 
low as 4 km/s, and material heated up to ~1200 K (petro¬ 
logic type 6) is expected for most impacts in the first 100 
Myrs (Davison et al., 2013). 

Numerical simulations of crater formation in large sub- 
catastrophic impacts also suggest that heated impactor and 
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Fig. 5.— The probability that a given mass fraction of a meteorite parent body will be heated to a range of temperatures, for (a) those 
asteroids that survive the first 100 Myrs of Solar System history without being disrupted and (b) those bodies that are disrupted. Modified 
from Davison et al. (2013) 
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target material, and exhumed hot material from the deep 
mantle, is rapidly buried beneath an insulating lens of 
porous breccia. The structure of these lobate/hemispheric 
deposits, for the specific case of the global scale impacts 
Rheasilvia and Vennenia on asteroid Vesta, is studied by 
Jutzi and Asphaug (2011) and Jutzi et al. 2013. Buried 
impact-heated or impact-exhumed material cools slowly, 
consistent with the wide range of metallographic cooling 
rates inferred for chondritic materials (Davison et al., 2012). 

While the volumetric significance of impact-heated ma¬ 
terial appears to be consistent with the limited constraints 
on meteorite thermal metamorphism, compared to Al de¬ 
cay it appears that localized impact heating is an impor¬ 
tant secondary effect, even when considering the cumula¬ 
tive effect of multiple impacts over the lifetime of a par¬ 
ent body. There is of course a bias, that the proximity of 
impact-heated material to the parent body’s surface makes 
it most accessible for excavation and meteorite delivery by 
subsequent collisions. 

3.3. Accelerated cooling by impact overturning 

If internal heating by short-lived radioactivity was domi¬ 
nant, then sub-catastrophic collisions would have played an 
additional important role in regulating thermal evolution. 
Ciesla et al. (2013) showed how a single large impact can 
have an enormous influence on the cooling of planetesimals 
heated internally by radioactivity. 

In their models (Figure 5), hot, thermally-softened mate¬ 
rial is brought up from the deep interior to flow out and over 
the surface of the planetesimal, turning the asteroid “in¬ 
side out.” The direct exposure to space of deeply exhumed 
molten silicates would lead to rapidly accelerated cooling 
and degassing; more generally, it is a considerable source 
of enthalpy for materials brought to the surface, something 
Asphaug et al. (2011) proposed as an origin of chondrules. 

In general, a much larger fraction of an internally-heated 
planetesimal experienced enhanced cooling than is heated 
by the impact, implying the nonintuitive result, that the ma¬ 
jor effect of sub-catastrophic collisions was to accelerate the 
loss of heat, rather than to deposit heat (Ciesla et al., 2013). 

The rapid cooling of iron meteorites is a longstanding 
problem (e.g. Wood 1964) that could perhaps be solved if 
melted core material buried deep inside a differentiated as¬ 
teroid could be exhumed by this mechanism, so that the iron 
might cool rapidly. But from many simulations (Love and 
Ahrens 1996, and thereafter) it appears that an iron core, 
once buried, is almost impossible to blast out by impact 
overturning. Furthermore it would sink back to the core, 
cooling there instead of at the surface. So the above pro¬ 
cess, of ‘breaking open the pie crust’, would affect mostly 
the thermodynamics of meteorites from the upper to mid¬ 
mantle of a melted planetesimal, not the core. 

It is possible, on the other hand, to achieve greatly ac¬ 
celerated cooling of iron planetesimals if we consider the 
exhumation of the projectile rather than the target. As 
noted above, planetesimal projectiles (but not targets) dur¬ 


ing SSCs can be ripped apart into beads of iron-rich ob¬ 
jects during hit and run collisions involving larger embryos 
(Asphaug et al. 2006; Asphaug 2010). The mechanism of 
mantle stripping of hit and run projectiles is discussed in 
more detail in the chapter by Scott et al. (2015), including 
its implications for the history and evolution of the asteroid 
parent body population. 

While mantle stripping by hit and run has been simu¬ 
lated dynamically in a number of studies, the thermal con¬ 
sequences of abruptly exposing interior materials to space 
awaits the kind of detailed study that Ciesla et al. (2013) 
have given to large scale cratering. Asphaug et al. (2011) 
consider exhumation in a more specific context to explain 
the cooling rates of chondrules, that are formed according 
to their model by inefficient accretion of incompletely dif¬ 
ferentiated melted planetesimals. For the general case, it 
is apparent that the accelerated-cooling scenarios defined 
above for mega-cratering would have comparable, more ac¬ 
centuated aspects if applied to the aftermath of hit and run 
collisions, that typically rip off half the mantle of a differ¬ 
entiated body without producing shocks (Asphaug 2010). 
Yang et al. (2007) invoke hit and run mantle stripping to 
explain the rapid cooling rates of the IVA iron meteorites, 
in a scenario that requires almost complete mantle strip¬ 
ping to leave behind only a thin rind of insulating silicates. 
This would require repeated hit and run collisions while the 
body was still molten - something that Asphaug and Reufer 
(2014) argue would have been probable for at least a few of 
the unaccreted bodies that became the progenitors of aster¬ 
oids. 

4. Rheology and disruption 

One goal of current numerical modeling is a unified sim¬ 
ulation of asteroid collision including fracture, crushing, di- 
latatancy, friction, and debris flow. The physical modeling 
borrows equally from dynamics, seismology, geomorphol¬ 
ogy, and impact cratering. Only recenty has it been possible 
to look at the additional effects of friction, porosity and co¬ 
hesion, at least on a global scale, so this represents a frontier 
of computational geophysics. 

Jutzi (2014) used SPH to model colliding bodies in an 
exploration of the threshold for catastrophic disruption. Be¬ 
low we show a subset of those calculations: colliding aster¬ 
oids that are the same size in each simulation {Rt — 100 
km, Rp = 27 km), and the same mass (initial bulk density 
p = 1.3 g/cm^), and the same impact parameter (vimp = 3 
km/s, 6 = 45°). What is different between the four simula¬ 
tions in Figure 6 are the material properties: 

1. no friction, no crushing, no cohesion 

2. friction, no crushing, no cohesion 

3. friction and crushing, no cohesion 

4. friction, crushing and cohesion 

For the two cases with crushing, a crush curve with pa¬ 
rameters for pumice is applied (Jutzi et al. 2009). For the 
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Fig. 6.— Post-impact temperature distribution within planetesimals internally heated by radioactive decay of ^®A1 for a single impact 
that occurs at different times during thermal evolution. At early and late times, the impact occurs on a relatively cool planetesimal, 
resulting in a recognizable crater with a localized, impact-heated region at its center. In contrast, impacts that occur during the period 
of high internal temperatures (~2-20 Ma for this size planetesimal) result in uplift and overturn of deeply-derived, hot and thermally- 
softened material. From Ciesla et al. (2013) 




Fig. 7.— Cross-sections of four SPH simulations of collisions between targets Rt = 100 km and projectiles Rp = 27 km, at relative 
velocity 3 km/s and 45° impact angle. Four different rheologies are investigated. The degree of disruption, the size of the largest remnant 
(Mir/Mtot), and the increase in density strongly depend on the target properties. From Jutzi (2014). 
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case with cohesion, an internal friction = 1.5 and cohe¬ 
sion Fq = 100 MPa are applied, which is reduced to Hd = 
0.8 and = 3.5 GPa following damage. The first case in 
Figure 6 is a fluid, the second is a fluid resisted by friction, 
the third is a crushable fluid resisted by friction, and the last 
is a cohesive porous body like pumice. 

Based on a suite of such models, Jutzi (2014) found the 
catastrophic disruption energy Q’^ to be ~ 5—10 times 
higher when friction is included. The disruption thresh¬ 
old further increases by a factor ~ 2 — 3 when the energy 
dissipation by compaction (pore crushing) is taken into ac¬ 
count. (Surprisingly, there is not much difference between 
the case with or without cohesion.) This implies that many 
events previously thought of as catastrophically disruptive, 
are more likely global scale impacts that leave behind a 
completely transformed asteroid, but do not disperse it. 

These simulations are consistent with previous models, 
that cohesion and tensile strength are negligible in poten¬ 
tially disruptive impacts involving bodies larger than a few 
100 m diameter. But it is a misnomer to speak of this as 
the ‘strength-gravity transition’, because gravity is not the 
controlling aspect for the response, either. For example, 
collisions involving partially melted silicate targets up to 
~ 1000 km diameter may be limited by viscosity (Asphaug 
et al. 2006). For solid asteroids and planetesimals up to 
100s of km diameter that have been beaten into rubble, fric¬ 
tion and compaction are expected to dominate, according 
to these calculations. Although gravity cannot be ignored, 
there may be no simple ‘gravity regime’ for asteroid im¬ 
pacts. 

5. Conclusions 

Global scale impacts on asteroids and their parent bodies 
are in many respects similar to the giant impacts that created 
the terrestrial planets, and there is great scientific value in 
making the connection. But at asteroid scales, we cannot 
ignore the internal stresses and their associated rheological 
responses - the breaking of rock, the flow of rubble, the 
opening of faults, the seismological vibration - in addition 
to the self-gravitational complexities. 

Gravity is subtle, and so are the other forces, and approx¬ 
imations or even small numerical errors can overwhelm the 
accurate calculation of an asteroid impact. In recent years 
computer modeling has begun to catch up with the physics, 
and progress on this front defines much of the modern sci¬ 
ence of global scale impacts. In the coming decade, we 
can look forward to ever more sophisticated simulations; 
for now it is too early to say whether the models are suffi¬ 
cient to be predictive. 

Meanwhile, every new mission to an asteroid - and in 
particular, every new mission to a different size or class of 
asteroid - has revealed a solar system process not previ¬ 
ously envisioned, and has turned our previous understand¬ 
ing of small bodies on its head. Flyby missions are of great 
value, as are the increasingly impressive radar shape models 
(Ostro et al. 1990) and adaptive optics models (see chapter 


by Marchis et al. 2015) obtained so far for dozens of NBAs 
(see chapter by Benner et al. 2015). While lacking the fi¬ 
delity of a spacecraft flyby, radar images are sufficiently 
detailed to play a fundamental role in resolving the geo¬ 
physical debates that have been summarized in this chapter. 

Rendezvous missions provide a much more thorough ex¬ 
amination of a smaller subset - three asteroids so far, with 
three more on the near horizon; the sample return missions 
Hayabusa 2 and OSIRIS-REx to small NEOs, and the Dawn 
mission continuing to Ceres. Hayabusa 2 will conduct the 
first impact experiment on an asteroid, and although it is 
almost certainly too small of a scale to result in any of the 
global scale phenomena described above, it will no doubt 
shed light on some of the wide-open questions concerning 
microgravity impact physics at several meters scale. 

Efforts are ongoing to study the effect of a small space¬ 
craft colliding into an asteroid using an LCROSS-like mis¬ 
sion approach (Colaprete et al. 2010), which depending 
on the asteroid could cause global accelerations greater 
than local gravity, rearranging the global topography. Such 
a mission would elucidate the many factors that dictate 
whether a cratering event has global consequences or is lo¬ 
cal to one hemisphere - an understanding which in turn will 
allow us to put the more readily attainable cratering obser¬ 
vations from radar and flyby missions into context. 

And the last front is meteorites, including the samples 
that are collected by missions. Meteorites have a way 
of presenting us with seemingly contradictory information, 
leading to a spectrum of contested hypotheses for their ori¬ 
gin. But there is a chance, through collisional and thermal 
modeling and geochemistry, to understand their most ba¬ 
sic physical properties and the processes they have experi¬ 
enced, from the formative collisions to the catastrophic dis¬ 
ruptions that have populated the Main Belt and near-Earth 
space. 
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